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Abstract 

In the problem of planetary formation one seeks a mechanism to 
gather small dust particles together into larger solid objects. Here 
we describe a scenario in which turbulence mediates this process by 
aggregating particles into anticy clonic regions. If, as our simulations 
suggest, anticyclonic vortices form as long-lived coherent structures, 
the process becomes more powerful because such vortices trap par- 
ticles effectively. Even if the turbulence is decaying, following the 
upheaval that formed the disk, there is enough time to make the dust 
distribution quite lumpy. 

1 Introduction 



To rationalize the distribution and the motions of the planets and 
their satellites in the solar system, Kant 26 and Laplace 29 separately 
introduced forms of the nebular hypothesis, a conjecture that the solar 
system formed from a flattened gas cloud or disk. Their prescience 
has been confirmed: disks do occur in many places in the universe. 
Since much of the interesting behavior in disks has been ascribed to 
turbulence, a problem motivated by them is appropriate in a volume 
dedicated to R.H. Kraichnan. That problem is to understand how the 
planets formed in the primitive solar nebula. 
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Stars condense from an interstellar medium consisting mostly of 
gas with an admixture of solid particles called interstellar dust. Both 
observational and numerical studies suggest that, as the central star 
contracts, it leaves around it material that contains a good share of 
the initial angular momentum of the whole system. In this nebula, 
the centrifugal force balances the stellar gravity in the radial direction 
and a disk is formed. As the dust settles toward the midplane of this 
disk it is somehow accumulated into protoplanetary objects. 

The general picture has been reviewed by Lissauer 30 who has pro- 
vided an extensive bibliography. Other summaries exist (such as 
Ref.49) and there are a number of introductory expositions (for ex- 
ample Ref.ll). Two main possibilities have been considered. If the 
disk is not turbulent, solid dust particles settle into a very thin layer 
and then can accumulate into larger solid objects by gravitational 
instability 39 ' 24 . On the other hand, if the disk sustains turbulence, it 
has been thought that one must rely on a slower process of adhesion or 
growth by two-body collisions. We examine here an alternative pos- 
sibility in which turbulence in the early history of a protoplanetary 
disk leads to a rapid aggregation of particles by anomalous diffusion. 
Even if the turbulence is decaying, there is time for the particles to 
aggregate significantly. 

The idea that the formation of planets in the Kant-Laplace disks, 
was somehow mediated by the interaction of turbulence and rota- 
tion was adumbrated by von Weizsaker. 48 Subsequently, experiments 
showed how a field of vortices of opposite senses can produce inho- 
mogeneous particle distributions. 12,13 The experimental evidence that 
anticyclonic vortices (whose vorticity is in the sense counter to that of 
the prevailing global motions) can trap particles is now extensive. 45 ' 47 

An application of the mechanism of anomalous diffusion by anti- 
cyclonic vortices is to the formation of planets, 6 ' 44 as seen in several 
numerical studies of Langrangian orbits in prescribed vorticity fields. 36 
In this paper, we consider flows of the kind that would arise in a Kep- 
lerian disk and follow the aggregation process they produce. Our pur- 
pose is to show how this aggregation process can work despite some 
objections to the notion (to be discussed below) that protoplanetary 
disks can sustain turbulence. 4 ' 5 Without turbulence, it is unlikely that 
vortices on disks form at all. 

Our way around this impasse is to note that the initial aggregation 
of dust could take place in the decaying turbulence that was likely to 
have occured in the formation of the nebula. Once this has happened, 
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the usual sequel of gravitational clumping can be presumed to take 
over. We shall present simulations to show how anticyclonic vortices 
form from sufficiently large initial perturbations of Keplerian disks (see 
also Refs.f 0,2,41), just as in earlier simulations with mild shears, 32,21 
and strongly promote the aggregation of dust. The question we then 
investigate is whether the vortices have time before they die away to 
strongly aggregate the particles. 

In the next section, we write down the equations of motion to set 
the notation in order. We turn then to a brief summary of relevant 
aspects of disks, after which we exhibit results of simulations of a 
Keplerian flow. These show how long-lived anticyclonic vortices form 
from initial finite disturbances. Then, from a study of the Lagrangian 
orbits of dust in this flow, we show how an initially homogeneous dis- 
tribution of particles is redistributed into strong local accumulations 
of particles. 



2 Equations of Motion 

To concentrate on vortex dynamics, we restrict ourselves to incom- 
pressible flow with 

V-v = 0. (1) 
We also assume constant viscosity so that the equation of motion is 
Dv 1 

— = — Vp - V$ + yAv , (2) 
Dt p 

where $ = —GM/r is the gravitational potential of the central mass, 
the protosun. 

The key quantity here is the vorticity, for which the evolution 
equation is 

Dui 

■jjj; = (Ml ■ Vv + vNjJ . (3) 

In a disk, the large scale flow is nearly two-dimensional and we 
adopt this simplification here. Then we may express v in terms of 
a stream function ip such that, in Cartesian coordinates, the velocity 
field is given by 

®t ^t-\ (4) 
dy ' dx J 

The vorticity then has only one component, u, and it is related to the 
streamfunction by 

u = V 2 V> • (5) 
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Hence equation (]3|) may be written as 

^ + J{%b,to) = v&uj (6) 

where J(ip,u) = d x ipd y ui — d x ujd y ip. 

In the absence of dissipation (y = 0), equation (0) admits an 
infinite number of conserved quantities, two of which are quadratic 
invariants; these are the kinetic energy E = 1/2 / (Vip) 2 dxdy and 
the enstrophy Z = 1/2 J (V 2 xjj) 2 dxdy. The conservation of these two 
quantities in the inviscid limit is what makes two-dimensional turbu- 
lence so different from three-dimensional turbulence at large Reynolds 
numbers. 7,27 In two-dimensional turbulence there is a direct cascade 
of enstrophy from large to small scales and an inverse cascade from 
small to large scales of kinetic energy. 17,38,28 

For freely decaying, barotropic turbulence in shear-free environ- 
ments, intense long-lived vorticity concentrations are observed to form 
after an energy dissipation time. These coherent vortices are charac- 
terized by a broad distribution of size and circulation and they con- 
tain most of the energy and the enstrophy of the system. 33,34 The first 
question we take up here is whether the strong Keplerian shears of 
accretion disks prevent this well-documented behaviour. 



3 Astrophysical Disks 

The purpose of this section is to provide a general idea about the 
background for the fluid processes that we study in the next sections. 
This is called for in a controversial subject such as planet formation 
since many readers of this journal may not have previously encoun- 
tered some of the issues raised. On the other hand, the occurence of 
disks of many sizes throughout the universe is so well documented that 
our discussion can safely start from the assumption that one of these 
once surrounded the protosun. We shall not go into the processes by 
which disks are thought to form but consider only the simplest case 
of a disk that is axisymmetric in the large and symmetric about the 
midplane z = 0. In the radial direction, the principal support against 
the gravitational pull of the protosun is provided by centrifugal force 
while, in the z-direction, the support is by a pressure gradient. If the 
ratio of these two support forces is large, the disk is quite thin, as we 
shall presume here. 
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Let the rotational speed at distance r from the central protosun 
be v. The centrifugal balance condition is 

v 2 = 2GM/r , (7) 

where M is the mass of the protosun. This describes a Keplerian flow 
that is stable to small perturbations since it has no inflection points 
and the Rayleigh discriminant, 

A = r- 3 d{r 2 fl) 2 /dr , (8) 

is positive. In astronomy, where this stability condition is often met, it 
is custumary to introduce the epicyclic frequency k = a/A, whose role 
in the dynamics is not wholly different from the buoyancy frequency 
of stably stratified fluids. Since the disk is thin, we may identify r 
with the radius in cylindrical coordinates. In our two-dimensional 
simulations, we retain this meaning of r and the form —GM/r 2 for 
the gravitational acceleration. 

If the flow remains laminar, there will be a slow viscous transport 
of angular momentum outward that allows a weak inward spiralling 
of the material onto the protosun. The two-dimensional version of 
such motion has been studied in the context of drag reduction theory 8 
where it is thought that a slow inflow can inhibit boundary layer sep- 
aration and reduce the drag on airfoils. This process is too slow for 
astrophysical purposes when it is laminar, but its turbulent analogue 
has been invoked for the case of astrophysical disks, 31 which are called 
accretion disks because of this inflow. In the turbulent case, it is sug- 
gested, the dissipation rate could be high enough to render some disks 
quite bright. But can we assume that disks are turbulent? 

Until recently, many took it for granted that Keplerian disks were 
turbulent because the Reynolds number in a disk is literally astro- 
nomical. On account of (0) and the relatively small viscosity, it was 
often assumed 15 that "The successive rings of gas in the medium will 
therefore have motions relative to one another, and turbulence will 
ensue." But simulations by Balbus and collaborators, 3 among oth- 
ers, have not shown nonlinear instabilities in Keplerian disks at the 
respectable Reynolds numbers they achieved. 

Apparently, the circular Keplerian flow is different from plane Cou- 
ette flow, where nonlinear instability is known to occur, and indeed 
does so in the calculations run by Balbus et al. Nevertheless, their con- 
clusion is not universally accepted 23 and there is an issue of principle 
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to be decided here. One technique that may be decisive is the method 
of energy stability which may prove helpful in deciding whether there 
can be sustained turbulence in nonmagnetic Keplerian disks. 

On the other hand, for the typical astrophysical conditions, the 
matter in disks is ionized and is likely to be magnetized. A preva- 
lent magnetic field can destabilize a shear flow that would otherwise 
have been stable 16 ' 43 ' 18 and Balbus and Hawley 5 have introduced this 
notion into disk theory to great effect to justify the assumption that 
turbulence does occur in accretion disks. Once started, this magne- 
toturbulence may be expected to produce a dynamo and to be self- 
sustaining along with the turbulence. In this way, the study of disk 
structure can be said to have been reduced to a previouly unsolved 
problem. The literature on this problem is of course quite large, in 
keeping with its importance in astrophysics. 

By contrast, in the case of the disk that is supposed to have 
spawned the solar system, it is thought that the matter was too cool 
to be ionized and so the flux was not frozen in Ref.46. We cannot 
then call on magnetic effects to render the solar nebula turbulent. 
The problem of whether such a disk is turbulent is therefore open, 
and there are at this time no real grounds for invoking turbulence. 

However that issue is decided, it is likely that most would accept 
the assumption that the initial conditions in the disk were turbulent, 
and this is sufficient for our purposes, even if that initial turbulence 
were decaying. Indeed, our results indicate that the decay lasts many 
rotation periods, and that is ample time for the vortical aggregation 
process we describe here to make good headway. 

4 Simulations with Keplerian Shear 

If disks are turbulent, they might be expected to follow the pattern of 
other rotating turbulent fluids 9 ' 25 ' 33 ' 42 ' 32 ' 21 ' 22 and form coherent vor- 
tices that last many turnaround times. Even for turbulent disks, this 
possibility is not generally accepted because the shears are so strong in 
Keplerian flows. In this section we study this problem to see whether 
vortices may form and survive for long times in Keplerian shear given 
a finite initial disturbance. A positive result would imply that even 
if the protoplanetary disk were nonlinearly stable, it would have pro- 
duced vortices early in its lifetime since it is likely to have been formed 
in a turbulent state. In this section we find that for initial perturba- 
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tions above some threshold, long-lived vortices do form, though they 
eventually decay. Then, in the next section, we investigate whether 
the vortices have enough time in their decay periods to have a signifi- 
cant effect on the particle distribution in a disk and so contribute the 
formation of planets. 

To see whether vortices can form in a Keplerian shear, we per- 
form numerical simulations of an incompressible, two-dimensional flow 
driven with a shear like that derived from (7). Into this flow, we in- 
troduce sizable perturbations such as one might expect in an initially 
turbulent situation and follow the formation of vortices that indeed 
survive for many rotation periods of the disk. This approach is like 
that used previously 34 to simulate vortex formation in decaying tur- 
bulence and it permits us to study the effect of strong shear in this 
problem, which turns out to be that anticyclonic vortices are heavily 
favored. 

4.1 Technical details 

The basic quantity that we compute is the vorticity, which we split into 
two parts, a specified (regularized) Keplerian flow and a disturbance 
to the basic flow. Thus we write, 



where ( is the dynamical quantity of interest. 

The imposed background Keplerian velocity profile is, in regular- 
ized form, taken to be 



where K and r$ are adjustable parameters. For r > ro, a conventional 
Keplerian profile Kj^fr is recovered while for r <C ro, the velocity 
field reduces to a solid rotation Kr/r^ 2 . The behavior at small r 
is introduced to avoid the singularity at the origin associated with a 
point mass. This profile has a vorticity field, ujk = }:-§p( r VK), given 



(9) 




(10) 



by 



= ^ {l - expR^] + 3 (±f ex P[ - ( ^]} 



(11) 



For r <C r , ujk — SK/r^ and for r 3> r , ujk — K/2r 3 / 2 . 
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Our simulations are based on a standard pseudo-spectral code with 
2/3 dealiasing and resolution 512 x 512 grid points in a square box 
with periodic boundary conditions [— tt, tt] x [—it, tt]. In the simulations 
discussed below, we have used the value r$ = 0.123, which corresponds 
to 10 grid points in the region of rigid rotation. We also adopt K = 
2.07, which gives the mean kinetic energy of the Keplerian shear the 
value Ek = 0.5. With these choices, we ensure a Keplerian profile 
with UJK oc 1/r 3 / 2 for the underlying vorticity and velocity fields in 
the range 0.3 < r < 2 (see figure 1). For r < 0.3, the disk rotates 
rigidly and for r > 2 the pure Keplerian profile is altered by the 
periodicity of the domain and the interaction with disk images from 
the periodic boundary conditions. At r = vr/6, the angular velocity 
is 2lok — 5.46, which gives a typical rotation time Tk. — 1.15. The 
simulations have been run up to a total time T = 20. 

When we introduce the decomposition (||) into (||), we obtain an 
evolution equation for the disturbance vorticity, £: 

^L + J( ip ,() = vA( + F (12) 

F = -J(iI> k ,Q-J(<P,uk) (13) 
where <p and are the streamfunctions associated with £ and lok 
respectively. The quantity F represents the effect on the disturbance 
flow of the imposition of a background Keplerian shear, which satisfies 
J(^k,^k) = 0. We have assumed that the dissipation acts only on 
the disturbance field £ and we have taken v = 5 x 10~ 5 . 

As initial condition, £(r,0), we select a narrow-band random vor- 
ticity field with energy spectrum 34 

un 

v ; [(m/n)k + k] m+n v 1 

We have taken ko = 10, m = 30 and n = 5. The value of E$ is chosen 
to control the relative energies of the Keplerian and of the disturbance 
flow, as specified below. The Fourier phases are randomly distributed 
initially between and 2ir. In order to avoid unrealistic periodicity 
effects, the disturbance has been limited to a region of radius r max = 2, 
so that no perturbations occur near the edges of the simulation box. 



4.2 Formation of anticyclonic vortices 

In our simulations, when the initial energy of a disturbance is less 
than about 10 -3 of the energy in the Keplerian flow, the vorticity 
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fluctuations are sheared away and the disk returns to its unperturbed 
velocity profile. This behavior is consistent with the linear stability of 
Keplerian disks. When we begin a calculation with larger perturbation 
energies than this, an initially random vorticity field forms itself into 
well defined vortices, much as in flows with little or no shear, except 
that our vortices are all anticyclonic. Cyclonic vortices do not form 
and, even if they are expressly introduced initially, the shear quickly 
destroys them. 

Figure 2 shows the vorticity field for a simulation in which the 
initial disturbance energy was chosen to be 4 x W~ s Ek- The vorticity 
is shown at times t = (figure 2a), t = 8 (figure 2b) and i = 15 
(figure 2c). The presence of anticyclonic vortices is clearly visible. 
Because they last many of their own turnover times, such vortices are 
an example of coherent structures. 

The initial sizes of the developing anticyclonic vortices are small, 
but they merge and grow until they apparently reach a limiting size. 
The limiting size of the vortices increases with the distance to the 
center of the disk roughly like r 1,5 out to r ~ 1.5, though small vortices 
penetrate into the central regions of the disk. In this system, formation 
and persistence of coherent vortices is possible only when the nonlinear 



term J((p,Q in eq. (12) is larger than dissipation z^A£, large-scale 
shear J(ipK, C) an d differential rotation J((f, ^k)- For Keplerian disks, 
the range of scales where vortices can persist is thus limited from 
below by dissipation and from above by shear and differential rotation. 
Whenever the two limiting scales become too close, vortices cannot 
survive. 

For thinking of the solar nebula, we identify a length scale such 
that the border of the disk, at r max = 2, corresponds to 7AU. (An 
astronomical unit is the mean distance of the earth from the sun.) 
With this convention, the biggest vortices of our simulation are located 
at typically 5AU, at the distance of Jupiter from the center, while the 
smallest vortices lie at roughly the earth's distance from the sun. 

Once the length scale has been chosen, the time scale is determined 
by Kepler's law of rotation. First we recall that the rotation period 
is Tk — 1.15 at r = tt/6 in adimensional units. In physical units, this 
distance corresponds to about 1.83^4f7. At this distance, the (physical) 
revolution time is T£ ~ 2.48 yr. Physical time is therefore related to 
time in adimensional units by 

t*(year) = 2.16 t . (15) 
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To give the value of the viscosity coefficient a physical meaning, 
we recall that turbulent viscosity in disks is usually expressed in terms 
of the dimensionless a-parameter 40 by 

v = aH 2 uj* K (16) 

where H is the thickness of the disk and ui* K is the Keplerian rotation 
in the region of interest (in physical units). In the standard model 
of the solar nebula, a is in the range 1CP 3 — 10~ 2 in the region of 
Jupiter, assuming H ~ 0.3Af\ This determines a physical eddy vis- 
cosity y~5x 1CT 5 — 5 x 10~ 4 (AU) 2 /year. In adimensional units, it 
corresponds to 8 x 1CT 6 — 8 x 10~ 5 . The viscosity used in the previous 
simulations, i/ = 5x 10 -5 , lies in this range, indicating that formation 
of coherent vortices is indeed possible when the parameters take the 
values commonly estimated for the solar nebula. 

5 The Aggregation Process 

Anticyclonic vortices trap small particles. 4 ' 6 ' 44 ' 36 In this section, we 
study the effectiveness of this process when an explicit solution of the 
equations of motion is used to provide the vortices. In particular, we 
ascertain that the slowly decaying vortices in Keplerian turbulence 
have time to produce significant particle aggregation and so initiate 
planetary formation before they die away. 

5.1 The Particle Equation 

We next follow the Langrangian motion of dust particles in the flow 
computed in the previous section. This motion is described in a rigidly 
rotating frame with an angular velocity that is approximately the 
mean of the Keplerian rotation over the disk, namely f2 = 0.617. 

The study of particle motion in a fluid flow is surprisingly complicated, 35 '' 
especially for fluffy, extended flakes. But here, we consider the sim- 
plest case of dust particles whose sizes are less than the mean free 
path in the ambient gas and whose typical density, p^ is much greater 
than p g , the density of the gas around it. For this case, the equation 
of motion of a dust particle with position r is 

;F = - MA *-H*- T(r ' t) ) + ( n — ) r (17) 
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where v(r, t) is the fluid velocity in the rotating frame and 7 is a fric- 
tion coefficient. In addition to the gravitational force on the particle, 
this equation contains the Coriolis force, the centrifugal force and a 
phenomenological drag force produced by the ambient medium. 

For the rarefied conditions of the protoplanetary nebula, 7 oc p g /a 
where a is the particle radius; this corresponds to the so-called Epstein 
regime. 30 Our simulations of the particle motion show that the value 
of 7 is important in the aggregation process. For large 7, the particles 
are rather like tracers and for small 7 the concentration process is 
more rapid. 

5.2 Trapped particles 

We start our simulation of the Lagrangian trajectories with a uniform 
distribution of particles moving initially at the local fluid velocity. 
This leads to the results shown in figure 3 for 7 = 10. The successive 
panels display the particle distributions at times t = 2,8, 15 and 20. 
To bring out the effect of vorticity, £, in figure 4 we provide a plot of 
the time evolution of the average density of particles in three different 
regions, namely where C < —Co (full line), —Co < C < Co (dotted line) 
and C > Co (dashed line) with Co as an estimate of the r.m.s. initial 
vorticity. The plotted densities are divided by the mean density, so 
that the initial distribution is everywhere unity. We observe that the 
particles are expelled from the cyclonic regions and concentrated in- 
side the anticyclonic vortices. Simulations with other values of the 
drag parameter 7 have confirmed this general picture; the particle 
concentration process is even more rapid for larger particles associ- 
ated with smaller values of 7. A full study of the role of the various 
parameters appearing in eq.(|l7]) is reported in Ref.44 for a kinematic 
flow field. 

These results are qualitatively consistent with earlier empirical 
work 47 ' 12 and with theoretical studies in which the particle paths were 
followed in an a priori prescribed flow containing vortices 6,44,14 . The 
quantitative result added here is that decaying vortices in a real flow 
derived from the fluid equations have adequate time to aggregate par- 
ticles. 

Thus, by marching the particle motion along with the simulations 
that produce the vortices we have been able to see how turbulent 
fluctuations initially present in the disk draw particles into regions of 
negative vorticity. Already in the early evolution of the solar nebula 
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we observe the development of pockets of dust which fill the disk in- 
homogeneously even as the vortices form. This separation of the dust 
into different regions is caused by the cyclone /anticyclone asymmetry, 
and it is so rapid that at t = 2 (a matter of years) the density of the 
dust particles in the region C < — Co has already been increased by a 
factor 2. When the vortices do form, they are already surrounded by a 
large number of dust particles and so the rate of capture is enhanced. 

6 Discussion 

Rotating cosmic bodies whose atmospheres are turbulent and can be 
seen with good resolution display coherent structures — vortices in 
the atmospheres of planets and magnetic flux tubes in the solar atmo- 
sphere. However, many accretion disk theorists, feel that this obser- 
vation is not a good indicator of what happens on disks because they 
feel that the strong Keplerian shears in accretion disks will destroy 
vortices. There is practically no direct observational evidence on this 
point, but indirect evidence does suggest that coherent structures form 
in the disks around quasars. 1 Our fluid dynamical simulations, (see 
also Ref.10) show that Keplerian shears feed anticylonic vortices and 
destroy only the cyclonic ones. For vortices to form however strong 
perturbations such as might be provided by turbulence are needed. 

It is believed that the progenitor of the solar system was so cool 
that is was not ionized and so magnetic flux would not have been 
frozen in. Then the powerful magnetic instabilities that render hot 
disks turbulent would not operate and there is no well established 
mechanism for sustaining turbulence. Nevertheless, we have argued 
here that one may call upon a residual turbulence from the formation 
throes of the solar nebula. The vortices that would have been formed 
in those initial stages do seem to have time to act significantly on 
the particle distribution before decaying away. It is even possible 
that the process might have involved some of the primordial fields 
in the collapsing interstellar matter. Such details aside, it does seem 
reasonable to conclude that, once there is some turbulence, anticylonic 
vortices form. This in itself is an interesting fluid dynamical process 
that is likely to take place on conventional accretion disks. The version 
of the problem that we have looked at has features in common with 
two-dimensional turbulence, which R.H. Kraichnan has done much to 
explicate. 27 ' 28 



12 



Of course, the situation is more complicated in the real astrophys- 
ical contexts than the version we have presented here and a more re- 
alistic version of the theory might be based on shallow layer theory. 10 
Even if if the fluid is compressible, as it is in astrophysical disks, 
the thin layer approximation is very close to the standard shallow 
water equations, 37 which also give rise to vortex formation if strong 
disturbances are introduced. 19 ' 20 Whether vortices form in disks or 
not, large vortical disturbances can by themselves redistribute parti- 
cles. This may in itself be enough to initiate planet formation, but a 
complete theory would require a more detailed knowledge of how the 
nebula formed. 

In more general terms, the (perhaps) surprising ability of turbu- 
lence to render an initially homogeneous distribution of ambient par- 
ticles very inhomogeneous is likely to have consequences in several 
fields and deserves further study in the wake of the recent theoretical 
efforts in this direction. 47 ' 14 ' 36 This problem might also profit from a 
statistical treatment analogous to what has been done in the problem 
of the advection of a passive scalar, another subject to which R. H. 
Kraichnan has made a significant contribution. 
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Figure Captions 

Figure 1. The background velocity profile used in the simulations. 
The conventional Keplerian profile is plotted as a dashed line and the 
regularized profile is in full line. 

Figure 2a. Initial vorticity field of the perturbed two-dimensional 

Keplerian disk (t = 0). 

Figure 2b. Vorticity field at t = 8. 

Figure 2c. Vorticity field at t = 15. 

Figure 3a. Distribution of the particles in the disk at time t = 2 for 
7 = 10. 

Figure 3b. Distribution of the particles at t = 8. 

Figure 3c. Distribution of the particles at t = 15. 

Figure 3d. Distribution of the particles at t = 20. 

Figure 4. Evolution of the average density of particles in different 

regions of the disk characterized by specific values of the vorticity: 

C < -0.5 (full line), -0.5 < ( < 0.5 (dotted line) and C > 0.5 (dashed 

line) . 
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